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Abstract 



A study of an algorithm method capable to reveal anisotropic solutions 
of general scalar-tensor gravitation -including non-minimally couplings- is 
presented. It is found that it is possible to classify the behavior of the field of 
different scalar-tensor theories in equivalence classes, with the same classifier 
function that was obtained in Priedmann-Robertson- Walker models. 
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Scalar-tensor gravitation have proved to be a useful tool in the understanding of early 
universe models. The first and best known case of such theories is Brans-Dicke (BD) gravity 
|Jl|], in which there is a coupling function w{(f)) equal to a constant; (p being a dynamical 
field related with the previous gravitational constant. More general theories with other 
couplings have also been studied 0. The interest on these theories have been recently 
rekindled by inflationary scenarios and fundamental theories which seek to incorporate 
gravity with other forces of nature. Particularly, in string theories, a dilaton field coupled 
to curvature appears in the low energy effective action [^]. When scalar-tensor gravitation 
is concerned, one is interested also in the cosmological models it leads. Observational con- 
straints, mainly coming from the weak field tests |0 and nucleosynthesis 0, put several 
bounds upon the couplings. In any case, in order to evaluate the cosmological scenario and 
to test the predictable force of any theory, it is desirable to have exact analytical solutions 
of the field equations. But it was not until a few years ago, that Barrow [0, Barrow and Mi- 
moso IP and Mimoso and Wands ||^ derived algebraic-numerical methods that allow exact 
Friedmann-Robertson- Walker (FRW) solutions to be found in models with matter content 
in the form of a barotropic fluid for any kind of coupling w{(j)). 

However, scalar-tensor theories can be formulated in two different ways depending on 
the choice of the basic action. These two possible choices are the BD one, in which there is 
an arbitrary function in the kinetic term of the scalar, and the one which admit an arbitrary 
function multiplying the curvature while mantaining a common kinetic term, for the theories 
known as non-minimally coupled (NMC). Via a field redefinition one can stablish the equiv- 
alence between these choices in the most of the cases, but that is not what happens when 
the functions involved are not analytically invertible as, for instance, in the hyperextended 



inflationary scenario [|T0[. Very recently, we have presented a study on the full lagrangian 



density for the field, which involves, in the more general case, two free functions [^. This 
lagrangian reads 

L = 167rL^-^0,^< + G(0)-ii? (1) 



where, as usual, Lm refers to the matter lagrangian density and R is the curvature scalar. 
Each possible choice of the action may be reproduced by a convenient selection of G and/or 
w. We have called hyperextended scalar-tensor gravitation (HSTG) to the theories of gravity 
this lagrangian leads, because of the similarity with the inflationary model. Similar algorithm 
methods of massless scalar fields, developed by Mimoso and Wands in were applied to 
this general approach. Those methods allowed to compute exact analytical solutions for 
FRW models in vacuum or with matter content consisting in radiation or stiff fluids for any 
choice of G and w simultaneously lTl|]. That includes the cases of NMC, where solutions are 
scarce . As spinoff of that research, we find that it was possible to define an equivalence 
classes behavior of scalar-tensor gravitation in which the scalar field itself is a class variable. 
Specifically, for different lagrangian densities -a set (G, w)- that have an equal function of 
given by 

we find the same solution for the field as a function of the conformal time Q |]TT|. Then, 
inside each class, the scale factor could be obtained inmediately as = GX where X was 
a class-dependent function of time. This approach was also used in the search of slow-roll 
solutions for non-minimally coupled theories |[T3[| . 

Anisotropic homogeneous cosmological models are being intensively studied since quite 
a long time [jl4[. In particular, within scalar-tensor gravitation, the analysis of anisotropic 
cosmologies could reveal different behavior when compared with Einstein General Relativity 
near the singularity or in the inflationary epoch |]16[. Processes of isotropization of 



Brans-Dicke Bianchi-type solutions are also of current interest [|T^]. 

The aim of this work is to explore, within the general lagrangian density of HSTG, the 
case of some anisotropic cosmological models. In particular, we are interested in examine 
if is it possible to extend previously derived results for anisotropic models in BD gravity 



The conformal time is defined as dr] = adt where a is the scale factor 
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to this general approach; allowing, for instance, to compute also anisotropic solutions for 
non-minimally coupled theories. In addition, we want to observe if anisotropy do or do 
not break the classification scheme with a as the classifier function. Throughout we draw 
heavily on the results on massless fields in anisotropic universes [|T8| and in the properties 
of the lagrangian (|I|) when considered in the Einstein frame. 



The field equations of HSTG are Q 



R^iv - -jp^vR = G{(j)) 



(3) 



1 doj uj ^ 



dG~ 



+ ^ 00 + 3G ^□(G-i) - G ^8nT = (4) 



It is very important to remark that the usual relation T^'^.j^ = stablishing the conservation 
laws (in the meaning of GR) of the matter fields holds true. Note that when G = l/4>, these 
equations reduce to the common BD ones. The Einstein frame was introduced by Dicke |]TU 



when working in Brans-Dicke gravity, by defining a conformal transformation in the form 



Qab — G. 



0(P9ab 



(5) 



where Go is an arbitrary constant which becomes the gravitational constant in the trans- 
formed frame. In a similar fashion, we introduce 



9ab = GoG{(t)) gab 



(6) 



Using the relation between the curvature scalars of the common -Jordan- frame and the 
transformed -Einstein- frame, given for all conformal transformations by Synge we get 
for the action 



1 
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1 5 ~ab(P,a4>,b 3 , r 



(7) 



Equation (Q) is obtained after the assumption that (f> depends only on time. 
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where we have defined Lm = Lm/ {GoG{(f)) ^)^. It is now possible to define a new scalar field 
^ by 



16nG, 



-a- 



(8) 



such that, 



S 



EF 



IGvr 



1 ~ 



1 



—R+lQn{Lm--^,at") 



(9) 



Thus, we recover the Einstein action with a stress energy tensor given by the sum of two 
contributions, the matter and the scalar field ones. It is worth noting that the scalar field 
is proportional to the variable Y which was introduced to solve the problem of FRW 
models in The recovery of the Einstein action -and the field equations derived from 
it- does, however, have a cost. That is, a non-independently conserved stress-energy tensor 
for matter, which in the Einstein frame behave in agreement with 

' 6 dG-^\ 



~ 1 IGttGoI ~ , / 



(10) 



So, in any case in which T ^ 0, only the total stress-energy tensor will be conserved, i.e. 
(Tab + T^y^ = 0; where stands for the stress-energy tensor related with the field, which 
is given by 



1p,a1p,b - -1p,ci''''gab ) • 



(11) 



Equation (|TT]) and cosmological assumptions for the field show that it behave as a stiff fluid 
with density and pressure given by: 

2 



1 (bp 
2 



(12) 



Due to the exact reproduction of the the Einstein fleld equations, the common results of 
general relativity will hold good. We shall take into account matter given by a barotropic 
fluid, particularly in the cases of stiff fluids or radiation or both. The case of a dust fluid plus 



radiation in scalar-tensor theories is analyzed in ||2l| while some cases of imperfect fluids in 



2^ . Following Raychaudhuri p3|, and considering those models in which the velocity of 



matter is parallel to the unit normal to the spatial hypersurfaces (a geodesic time-like vector 
it is possible to write the Einstein field equations -the HSTG equations in the Einstein 
frame- in the form of a constraint 

^2 = 247rGo(p + p'^) + 3a2-^ (13) 

plus the Raychudhuri equation. 

Hf) 1 

_ + -e' = -AnGoip + + 3(P + /)) (14) 

where we have introduced the expansion 9, the shear a and the curvature scalar of hyper- 
surface of homogeneity ^R\ all them, in the Einstein frame. The transformed quantities 
are 

2 

p = o p = o a' = Z ^ .^ dP = (GoG(4>)-') de (15) 

(GoG(0)-i)' ^ (GoG(0)-i)^ (GoG(0)-i) \ ^ ) ^ > 

we also introduce a volume factor V = {GQG{(f))^^ f^'^ V with V such that 6 = dV/Vdt. 

Having in hands both conservation laws, for matter in the Jordan frame and for matter 
plus field in the Einstein frame, it is possible to derive the corresponding energy densities 
and pressures. They are: 

Here, F is related with the presence of radiation and M with the presence of a stiff fiuid. 
Both, F and M, are positive constants. When a stiff fluid is present, it is possible to define 
a new field x, minimally coupled to the metric, such that its energy density be the sum of 
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the energy density of the scalar and the matter content. This field x can be related with 

0by iniin 



r- d(t) 

/a -f- (19) 



V 3 J V J 4MGoG(0)-i 

We see that, defining a new set of functions {G,w) by 

the effect of stiff fluid matter is equivalent to a new lagrangian with no matter content. 
That change is possible because the stiff fluid modify the dependence of x with (p, function 
that in vacuum models is solely accomplish by w and G. That was exactly what happened 
with BD models, where only existed w. Also here, the fleld equations become simplest with 
the use of the variable X = {GoG{(j))^^) = a} and the conformal time With this 



variable together with t/^ given above we can now respond to the two questions outlined in 
the begining of this work. To do so, let us treat in this formalism and as a flrst example, 
what happens with Bianchi I universes. 
The metric of Bianchi I models is 

ds^ = dt^ - ai{tfdx^ - a2{tfdy'^ - a-i{tfdz^ (21) 

Here, the expansion is given in terms of an averaged scale factor = 010203 = 1^ in the 
form: 9 = 3da/adt. The spatial curvature is null and the metric reduces to the flat FRW 
one in the case in which The general relativistic result for the shear holds 

in the Einstein frame; i.e. = 3S^/4a^, where S is a constant and d^ = V. Using the 
expressions for the energy densities, it is possible to obtain the constraint equation in the 
variable X. For matter given in the form of non-interacting stiff and radiation fluids it is 

= ^2 + ^2 + ATX (22) 

It has exactly the same form that the equation obtained in and so, it admits the same 
solution (equation (105) of that paper). Equation (^) is in fact a general relativistic result. 



valid here because of the Einstein frame [El]. It is important to stress that, although the 



equation and its solution are the same, the meaning of the variables are different and that 
it is now allowed the study of the more general case of lagrangian (|Ij). The solution of ( |22D 
shows a shear-dominated evolution at early times and a radiation-dominated evolution when 
the averaged scale factor tends to infinity. Note that these results do not depend neither on 
the particular form of G nor of w. The specification of X, the shear and each one of the 
scale factors of the metric -which can be obtained from the general relativistic results in the 
Einstein frame- describe the full evolution of the system. To go back to the Jordan frame, 
we have to obtain ^(r/). To do so, we have to invert (|19D and get and afterwards, use 
our knowledge of X{ri) to get xiv)- Both operations yields finally, (t){ri). From equation 



(pjp we see that to equal functional form of a^ac-, equal solution for the field is obtained. 
Indeed, it means that if a solution for in a particular BD gravity with w = w^ac is known, 
it is also a solution for the set of theories given by a^ac = {'^uiyac + 3)/3. Thus, it is possible 
to define a^ac as a classifier of equivalence sets. Within each one of these sets, the definition 
of X allows us to obtain the particular behavior of the averaged scale factor. The functional 
form of the metric scale factors does not depend on w or G in the Einstein frame, while 
it does in the Jordan frame. Is in this frame in which a particular dependence of G on 
discrime among different behaviors inside a equivalence set. As was the case in BD gravity, 
without the specification of the functions involved in the lagrangian, some conclusions may 
arise. In presence of a stiff fluid, a bounce will occur when 



which requires 

{A' + J:')a = {A'-AMGoG-')\^^^j . (24) 

It may be also seen that, if G^^ vanishes faster than X^, an anisotropic initial singularity 
in the Einstein frame becomes isotropic in the Jordan frame. Note that the condition for 
the bounce is incompatible with a vanishing G^^ and X. The possibility of having a finite 



expansion in this situation is analized for BD gravity in 



The above analysis implies that the algorithm method developed by Mimoso and Wands 
for Brans-Dicke theory is capable to deal also with Bianchi I models in general hyperextended 
theories and that the specific choice of a convenient w allows us to study also all the cases of 
non- minimally coupled gravitation which still retains the positivity of a, as was the case in 
BD gravity imposing, accordingly, values of w bigger than -3/2. The classification scheme 
of FRW models is not broken by anisotropy. 

As a matter of fact, one can -at this stage- observe that the same applies to others 
models, such as Bianchi V or III. Those models are studied in |^ for General Relativity 



and in |18[ for BD gravity, and the only thing necessary to answer the former questions in 
this same way is translate that analysis to the variables X, ip and equations (p!9D and (pO[) 
of this work. In the sake of conciseness we do not do so here. 
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